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Abstract 

' We present a new Skyrme energy density functional (EDF) named 
SAMi This interaction has been accurately calibrated to repro- 

' duce properties of doubly-magic nuclei and infinite nuclear matter. 

, The novelties introduced in the model and fitting protocol of SAMi 
are crucial for a better description of the Gamow- Teller Resonance 
(GTR). Those are, on one side, the two-component spin-orbit potcn- 

1 tial needed for describing different proton high-angular momentum 
spin-orbit splitings and, on the other side, the careful description of 

1 the empirical hierarchy and positive values found in previous analysis 

, of the spin (Go) and spin-isospin {G' ) Landau-Migdal parameters: 
< Go < Gq, a feature that many of available Skyrme forces fail 
to reproduce. When employed within the self-consistent Hartree- 
Fock plus Random Phase Approximation, SAMi produces results on 
ground and excited state nuclear properties that are in good agree- 

- ment with experimental findings. This is true not only for the GTR, 

, but also for the Spin Dipolc Resonance (SDR) and the Isobaric Ana- 
log Resonance (IAR) as well as for the non charge-exchange Isoscalar 
Giant Monopole (ISGMR) and Isovector Giant Dipole (IVGDR) and 

' Quadrupole Resonances (IVGQR). 

1. Introduction 

' The Skyrme Hartree-Fock (HF) plus Random Phase Ap- 
■ proximation (RPA) approach is one of the successful tech- 
niques for the study of the ground state and excited state 
properties of nuclei [SJ [3] . The Skyrme ansatz is also used 
in more elaborated theoretical frameworks that include 
higher-order nuclear correlations [2 13] . 

The Skyrme HF+RPA approach enables an effective de- 
scription of the nuclear many-body problem in terms of a 
local energy density functional. Specific drawbacks and 
problems exist [B] in this type of functionals. Therefore, 
we need to understand the origin of such deficiencies and 
eventually solve them. One of these problems is the ac- 
curate determination of the spin-isospin properties of that 
functional. If achieved, it should lead to accurate predic- 
tions of the properties of Gamow Teller Resonance (GTR) 
[7J, namely, the main focus of this work for its relevance 
in electron-capture in the core-collapse of supernova? [HUH]; 
neutrino-induced nucleosynthesis [TU| ; the study of double- 
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j3 decay and, if ever observed, for determining the 
neutrino mass in ncutrinolcss double-/? decay. GT matrix 
elements are also very useful for the calibration of detec- 
tors aiming to measure electron- neutrinos |12j . 

The sum of all possible transitions from the ground state 
|0) to any possible excited state \v) give the linear response 
or GT strength function, 

R GT ±(E) = E IH|0gtI|O>| 2 <5(£ - E v ), (1) 

V 

where Oq T is the GT operator 

A 

OtT = E <T W T ±W> ( 2 ) 

i=l 

and A is the mass number and <x and r are the spin and 
isospin Pauli matrices, respectively. The dominant tran- 
sitions will be those between spin-orbit partner levels as 
illustrated in the schematic picture in Fig.[TJfor the case of 
90 Zr. If this is true also for other nuclei, it is easy to real- 
ize that spin-orbit splittings between nucleon levels above 
- but still close to - the Fermi surface will play a relevant 
role. In this respect, it is well known in the field that most 
of the Skyrme interactions overestimate the experimental 
spin-orbit splittings in heavy nuclei |13j . 
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Fig. 1: Schematic picture of the most important single-particle 
transitions (o"t_) involved in the Gamow Teller Resonance of 
90 Zr. 

GTR measurements shows that the Ikeda Sum Rule 
(ISR) given by f[R GT - (E) ~R GT + (E)]dE = 3(N-Z) ex- 
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hausts only a 60-70% in the resonance region. To explain 
this quenching, it has been proposed that 2-particle 2-hole 
(2p-2h) correlations, or the coupling with the A— hole ex- 
citation, have to be taken into account. Recent experi- 
mental analysis in 90 Zr [TJ] seems to indicate that most of 
the quenching (around 20% of ISR) has to be attributed 
to 2p-2h coupling and that the role played by the A isobar 
is much less important. 

2. Skyrme Hamiltonian and Random Phase Ap- 
proximation 

The Skyrme interaction is a zero-range, velocity- 
dependent interaction that describe nucleons with space, 
spin and isospin variables. The Hamiltonian density is 
written as in Rcf. [15] . 

H = K + Ho + H 3 + H cS + Han + Hso + H sg + U Co ui (3) 

where K, = h 2 T/2m is the kinetic energy term where 
t = Y^t Pi j the zero range term, Hz the density de- 
pendent term, H c s the effective mass term, Hn n the finite 
range term, Hso the spin-orbit term, H ss is a term due 
to the tensor coupling with spin and gradient and Hcoui 
is the Coulomb term. The exchange part of the Coulomb 
term is calculated within the Slater approximation |16j . 
Specifically, 

n a = ^t Q [(2 + x Q )p 2 -(2x + l)(p 2 n + p 2 p )] 

n 3 = ^t 3 p a [(2 + x 3 )p 2 -(2x 3 + l)(p 2 l + p 2 p )] 

H eS = l[t 1 (2 + x 1 )+t 2 (2 + x 2 )} T p + 
o 

i [t 2 (2x 2 + 1) - ti(2x! + 1)] (r„p n + T p p p ) 

Hfin = ^ {[3^(2 + 2! )-i 2 (2 + x 2 )](Vp) 2 - 

[3* x (2a:i + 1) + t 2 (2x 2 + 1)] [(Vp„) 2 + (Vp p ) 2 ] } 
Hso = ^W J ■ Vp + i^(J ■„ Vp„ + J p ■ Vp p 

H Bg = - — (tiX 1 +t 2 x 2 )J 2 + —{t 1 -t 2 )(Jn 2 + Jp 2 ) 

The main differences with other Skyrme models avail- 
able in the literature are that we include a two parameter, 
Wo and Wq, spin-orbit potential J7]. For Wo = Wq one 
recovers the most standard form of this type of functional, 
for Wq = the spin-orbit potential mimicks relativistic 
calculations and for non- vanishing Wo ^ Wq, a mixed be- 
havior is found. And that we also include H sg known as 
central tensor term or J 2 term. 

Pairing correlations, important for the description of 
open shell nuclei, and deformation are not included in 
our calculations since we study closed shell, double-magic 
spherical nuclei. The center-of-mass correction adopted in 
our calculations uses the factor mAj (A — 1) instead of the 
bare nucleon mass m in the kinetic energy term K.. This 
prescription accounts for a large part of the center-of-mass 
correction and it is implemeted in most of Skyrme-HF cal- 
culations. 

The discrete RPA method we adopt in our calcula- 
tions is well-known in textbooks fTS] HH|- In our sen?- 
consistent approach, we build the residual interaction for 



the proton-proton, neutron-neutron and proton-neutron 
channels from the Skyrmc-HF energy density functional. 
Then we solve fully self-consistcntly the RPA equations by 
means of the matrix formulation. For further details see 
Rcf. [20] where, very recently, the code we have used for 
the calculations presented here has been published. 

3. Fitting procedure 

We present a non- relativistic functional of the Skyrme type 
PJ, named SAMi for Skyrme- Aizu-Milano. It is as accu- 
rate as previous Skyrme models in the description of nu- 
clear matter properties and of masses and charge radii of 
double-magic spherical nuclei. In addition, due to our fit- 
ting protocol, it improves the description of the GTR in 
medium and heavy mass nuclei with respect to previous 
models of the same type. 

For the minimization, we have performed a \ 2 test (see 
Ref. [JJ for further details) by means of a variable metric 
method included in the MINUIT package of Ref. [2JJ. We 
have chosen the following set of data and pseudo-data for 
our fit: 

(i) the binding energies of 40 < 48 Ca, 90 Zr, 132 Sn and 208 Pb 
and the charge radii of 40 > 48 Ca, 90 Zr and 208 Pb; 

(ii) the spin-orbit splittings of the lg and 2f proton levels 
in 90 Zr and 208 Pb; 

(iii) the fixed values for the Landau-Migdal parameters 
G = 0.15 and G' Q = 0.35; 

(iv) pseudo-data, corresponding to the variational calcula- 
tions of the energy per particle of uniform neutron 
matter at baryon density p between 0.07 fm~ 3 and 
0.4 fur 3 of Ref. [22]. 

The novelties of our protocol that guide the fit to a better 
description of the GTR are twofold and due to points (ii) 
and (iii). Actually, the impact on the excitation energy 
and strength of the GTR of (ii) and (iii) have been stud- 
ied in previous literature [23] [2H [25] . The importance of 
(ii) has been already discussed in the introduction. The 
hierarchy and values that we have taken for the spin and 
spin-isospin Landau-Migdal parameters in point (iii) have 
been empirically suggested by the study of Ref. [25] (find 
more details also in Ref. [TJ). 

On the other side, point (i) has allowed us to determine 
the saturation energy (eoo), density (poo) and incompress- 
ibility (Koo) of symmetric nuclear matter - constrained to 
be 240 ±20 MeV by an analysis of a large set of Skyrme in- 
teractions |27j - and point (iv) has been helpful in driving 
the magnitude (J) and slope (L) of the nuclear symmetry 
energy at nuclear saturation density towards reasonable 
values [H [2J [501131] 131 133] • In Tables 1 and 2 of Ref. [TJ 
the references for the used data and pseudo-data with the 
corresponding adopted errors, partial contributions to the 
X 2 , and the number of data points used in the fit can be 
found (Table 1); and the SAMi parameter set and some 
saturation properties with the estimated standard devia- 
tions (Table 2). 
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4. Results 

In this section we will show the results of SAMi as com- 
pared to other interactions and experimental data for the 
GTR in 48 Ca, 90 Zr and 208 Pb and the Spin Dipolc Res- 
onance (SDR) in 90 Zr and 208 Pb. Other important re- 
sults such as binding energies, charge radii and the com- 
parison of pure neutron and symmetric matter equations 
of states (EoSs) predicted by SAMi with state-of-the-art 
Brucckner-Hartree-Fock and the fitted variational calcu- 
lations of Ref. OH] can be found in Ref. [JJ . In addition to 
that, we have also tested that SAMi EoS is stable against 
spin and spin-isospin instabilities |25j up to a baryon den- 
sity of more than four times saturation density, well above 
the region important for the description of finite nuclei 
and enough for the study of uniform neutron-rich matter 
in neutron stars. 

The predictions of SAMi for important non charge- 
exchange excitations such as the Isoscalar Giant Monopole 
Resonance (ISGMR), Isovector Giant Dipole Resonance 
(IVGDR) and the Isovector Giant Quadrupole resonance 
(IVGQR) in 208 Pb are accurate. Specifically, the excita- 
tion (E x ) - or centroid (E c )- energies as predicted by our 
HF+RPA calculations (experimental values) for these res- 
onances are £ C (GMR) = 14.48 MeV (£ C (GMR) = 14.24± 
0.11 MeV [31]), -Bc(GDR) = 13.95 MeV (£ C (GDR) = 
13.25 ± 0.10 MeV [35]) and E x (GQR) = 23.1 MeV 
(^(GQR) = 23.0 ± 0.2 MeV [35]), respectively. The 
corresponding energy weighted sum rules agree well with 
available experimental data. 

The earliest attempt to give a quantitative description 
of the GTR was provided by the Skyrme SGII interac- 
tion [23j . Later on, an accurate functional for the pre- 
dictions of finite nuclei and charge- exchange resonances 
was proposed: namely SkO' [UJ. Relativistic mean- 
field and relativistic HF calculations of the GTR have 
also become available meanwhile [42] [43]. For this rea- 
son, we show in Fig. 03 the results for the Gamow- Teller 
strength distributions in 48 Ca (upper panel), 90 Zr (mid- 
dle panel) and 208 Pb (lower panel) as measured in the 
experiment [37] HU [35] [3jJ] HO] and as predicted the non- 
relativistic SkO' [41], SGII [23] and SAMi interactions ob- 
tained within our HF+RPA calculations. We also show 
the predictions of SLy5 [15j since it is of common use in 
the field and because our fitting protocol is similar to the 
one used to obtain this interaction and the predictions of 
PKOl [43] for the case of 208 Pb since it is based on a 
relativistic framework. 

In the upper panel of Fig. 03 we show the experimental 
data as well as the predictions of the SAMi, SGII and SkO' 
functionals for 48 Ca (the SLy5 result is not shown because 
RPA calculations produce instabilities, i.e., imaginary en- 
ergies). The excitation energy found in the experiment 
for the high-energy peak, -E° xp = 10.5 MeV, and that pre- 
dicted by SAMi, E^ AMi = 10.2 MeV, are in very good 
agreement. For the low-energy peak, the agreement of the 
excitation energy found in the experiment, E x xp = 3.0 
MeV, with that of SAMi, E* AMi = 2.0 MeV, is also good. 
The corresponding % of the ISR exhausted by the main 
peak between 5 and 17 MeV, is around 46% in the exper- 
iment and 71% in our theoretical calculations with SAMi 
(note that RPA does not include 2p-2h couplings). The 
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Fig. 2: GT strength distributions in 48 Ca (upper panel), 90 Zr 
(middle panel) and 208 Pb (lower panel) as measured in the 
experiment [37] [H EU \M E2 and predicted by SLy5 [T5], 
SkO' gj, SGII [23] and SAMi forces. For the case of 208 Pb we 
also show the predictions of PKOl [43] , 



overall description of the GTR in 48 Ca by SAMi improves 
the results obtained with SGII and SkO'. 

The prediction of the SAMi interaction in the case of 
90 Zr (middle panel of Fig. [2]) is even better than in the 
case of 48 Ca. The excitation energy and % of the ISR ex- 
hausted by the high- and low-energy peaks in the exper- 
imental data (SAMi) are, respectively, E° xp = 15.8 ± 0.5 
MeV and 57% (E^ AMi = 15.5 MeV and 70%) between 
12 and 30 MeV and E c x xp = 9.0 ± 0.5 MeV and 12% 
(#SAMi = 7 8 MeV and 27%) between 3 and 12 MeV. 
Again, the best overall picture is given by SAMi when 
compared with the other Skyrme interactions SGII, SkO' 
and SLy5. 

With unprecedented accuracy in HF+RPA calculations, 
the SAMi functional perfectly reproduces the excitation 
energy of the experimental GTR in 208 Pb [39] (lower panel 
of Fig. 02): E cxp = 19.2 ± 0.2 MeV and £f AMi = 19.3 
MeV. We also compare our results with the predictions of 
SGII, SLy5, SkO' and PKOl. None of them reproduces 
the excitation energy of the GTR in this nucleus within 
the experimental accuracy. 

The SAMi reliability for the calculation of other charge- 
exchange resonances such as the SDR or the Isobaric Ana- 
log Resonance in some doubly-magic spherical nuclei have 
been tested in Ref. [Tj. Here, we show the results for the 
SDR in 90 Zr and 208 Pb. The operator used for the RPA 
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Fig. 3: Spin Dipole strength distribution for Zr in the r_ 
channel as a function of the excitation energy E x measured in 
the experiment [S] and predicted by SAMi. Multipole decom- 
position is also shown. A Lorentzian smearing parameter equal 
to 2 MeV is used, 
calculations is, 

A 

EEriWrf [r L (f 4 )®<x«U (4) 
i=l M 

and, as it is shown in both figures, it connects single 
particle states differing by a total angular momentum: 
J" = CP, 1~ and 2~. The sum rule, 

J [R SD - (E) - R SB AE)]dE = ^(N(r 2 n ) - Z(r 2 p )), (5) 

is completely exhausted in our calculations, 99.99% in the 
case of 90 Zr and 100% in the case of 208 Pb. 

In Figs.[3]and|4]the SDR in 90 Zr and 208 Pb are compared 
with the experimental results for the total and multipole 
decomposition (J 71- = _ ,1~ and 2~) contributions. The 
overall agreement is noticeable. For the case of 90 Zr, the 
experimental value for the sum rule is 148±12 fm 2 [H] and 
the SAMi prediction is 150 fm 2 . The total and multipole 
decomposition (J 71 " = — , 1 _ and 2~) of the experimental 
data [ID] for the integral of i? SD - in the case of 208 Pb are 
1004t 2 ^ fm 2 , 107±7 fm 2 , 450^ fm 2 and U7±{ 6 5 fm 2 ; and 
for the SAMi predictions 1224 fm 2 , 158 fm 2 , 423 fm 2 and 
643 fm 2 , respectively. 

5. Summary and conclusions 

The new Skyrme interaction presented here and in more 
detail in Ref. [1] accounts for the most relevant quantities 
in order to improve the description of charge-exchange nu- 
clear resonances, i.e., the hierarchy and value of the spin 
and spin-isospin Landau-Migdal parameters and the pro- 
ton spin-orbit splittings of different high angular momen- 
tum single-particle levels close to the Fermi surface. As a 
proof, the GTR in 48 Ca, 90 Zr and 208 Pb and other charge- 
exchange resonances [T] are predicted with high accuracy 
by SAMi without compromising the description of other 
nuclear observables. 
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